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fT^ ^ In this paper a theorem is derived in order to provide a wide sufB- 

f^ ■ cient condition for an orthogonally transitive cylindrical spacetime to be 

^^D ' singularity- free. The applicability of the theorem is tested on examples 

(«_^ , provided by the literature that are known to have regular curvature in- 

O^ ■ variants. 
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I Introduction 

The issue of establishing whether a Lorentzian manifold is geodesically complete 
is not in principle a simple one since there is no Hopf-Rinow theorem that could 
settle the matter as it happens in the Riemannian case. One could think that 
regularity of the curvature invariants might be helpful, but there are known 
examples of spacetimes with regular invariants, such as Taub-NUT, |l||, that 
enclose geodesies that are not complete in their afHne parametrization due to a 
phenomenon called imprisoned incompleteness. 

Taking completeness of causal geodesies (^-completeness) as a definition of 
absence of singularities (no observer in free fall leaves the spacetime in a finite 
proper time), one can resort to many theorems in the literature (cfr. [n|, B and 
references therein) in order to determine whether a spacetime is singular. But 
on the contrary theorems that provide large families of nonsingular spacetimes 
are not very usual [|| , H and in principle the proof of geodesic completeness 
involves cumbersome calculations 0] . 

Instead of dealing with general Lorentzian manifolds, we shall approach or- 
thogonally transitive cylindrical spacetimes H since they have provided many 
examples of regular manifolds (cfr. H, m) in inhomogeneous cosmology H. 

Our aim will be the generalization of the theorem on diagonal orthogonally 
transitive cylindrical spacetimes in B to nondiagonal models and thereby com- 
prise all known nonsingular cylindrical perfect fluid spacetimes in the literature. 

First we shall show that the second-order system of geodesic equations can 
be reduced by the use of constants of motion to three first-order equations plus 
two quadratures. This fact will simplify the analysis of the prolongability of the 
geodesies and will enable us to write a sufficient condition for completeness of 
orthogonally transitive cylindrical spacetimes in a theorem. 

II Geodesic Equations 

We shall write the metric of an orthogonally transitive cylindrical spacetime 
in a chart using isotropic coordinates t, r for the susbspace orthogonal to the 
orbits of the isometry group and coordinates </>, z adapted to the commuting 
generators of the group of isometrics. The metric shall be determined by four 
functions, g, /, A, p, of the coordinates t, r, 

ds^ = e^ait.r) |_^^2 + dr2} + p-'it, r)e2/(*''-)d02 ^ 

+e-2/(*:-){rfz + A(r,t)d0}2. (1) 

and we shall assume that these functions are C^ in their range, 

- oo <t,z < oo, < r < oo, 0<(j)<2n. (2) 

The axis will be located where the norm of the axial Killing field vanishes, 
= A - .9(^,0 = p2(i,r)e2/(*^'-) +e-2/(*>'-) A\r,t), (3) 



which means that both A and p must vanish on the axis, since / is a smooth 
function. 

Since the choice of isotropic coordinates is not unique, we can take advantage 
of this freedom to have r = as the equation for the axis. In order to avoid 
conical singularities, the usual requirement llG] will be imposed in order to have 
a well defined axis, 

hm g(g'-^dA,gradA) ^ ^ 

r^O 4 A ^ ^ 

Denoting by a dot differentiation with respect to the afRne parameter, two 
of the four second-order geodesic equations, 

x' + T]^x^x'' = 0, (5) 

can be integrated taking into account that there are two first integrals of the 
geodesic motion associated with the generators of the isometries. These are the 
angular momentum around the axis, L, and the linear momentum along the 
axis, P, of a test particle of unit mass, 

L = e^ ■^(*^'-) p'^it, r) + e-2 /(*''■) A{t, r) {z + A{t, r) <^}, (6) 



P = e-2/(*'''){i + A(i,r)0}. (7) 

The affine parametrization is determined, up to an afhnity of the real line, 
by the prescription, 

S = e^sit.r) |p _ ^2| _ 1^ _ PA{t,r)yp-\t, r)e-2/(t,r) _ p2^2fit,r)^ (8) 

where S is one for timelike, zero for null and minus one for spacelike geodesies. 
Since we are dealing just with causal geodesies, for our purposes S will always be 
positive. After writing i, as functions of i and P, the second-order equations 
in t and r, 



i + (7t(i, r)i2 -I- 2 g,(i, r)if -I- gt (t, r)r2 _ p2g2 {/(t,r)-3(t,r)}^^(^^ ^) ^ 

~PAit,r)r l pt{t,r) PAtit,r) 

p2(i,r) \pit,r) ^^*^' >^ L-PAit,r 



^.,if(t,r)+,(t,r)y {L~PA{t,r)}^ f ptjt^r) ^ ^^^^^^^ ^ PA{t,r) ^^^ 



r + gr(t,r)P + 2gt{t,r)ir + gr{t,ry+P'e'if^'^^'>-i>^'^^»fAt,r)- 
^-2if(t,r)+,it,r)} {L-PAit,r)}^ f pr{t,r) , . .^ ^ , PAr{t,r) \ 

pHt,r) \pit,r) +^'-^*'"^ + L-PA(i,r)/' ^'"^ 

can be rearranged in a form that will be useful afterwards. 



{e29(t.Ot}- _ e-23(t,r)p^^^ ^)p^^^^ ^) ^ Q^ (^j) 



{e2ff(t,r)^|. ^p-2g(t,0i^(i^r)i^^(t^j,) ^ Q, (12) 



F{t, r) = e9(*''') \ S + P2e2/(t,r) + {^ _ p^^^^ ^)}2 ^^^ ^ (13) 

V p^i,^) 

that have the same estructure as in the diagonal case. 

If at least one of the constants L, P, 6 is different from zero, the system is 
equivalent to three equations of first order for future-pointing (past-pointing) 
geodesies, 

i^±e-'^s^'^''^F{t,r)cosh(,{t,r), (14) 

f = e-29(*''-)F(i, r)sinh e(t,r), (15) 

i(t, r) = -e-29(*''-) {±Ft(t, r) sinh^(i, r) + F,.(t, r) coshC(i, r)} , (16) 

parametrizing (jg) by a function ^(t, r). More explicitly, the last equation takes 
the form, 

a/(5 + A2p-2e-2/ + P2g2/ L V P / 

±smh^{Sgt+A^s^ht + P^e^fqt)}, (17) 

/i = g-/-lnp + ln|A|, A = L - PA, q = g + f, (18) 

that will be useful for deriving prolongability conditions for causal geodesies. 
The minus (plus) sign corresponds to past-pointing (future-pointing) geodesies. 
Note that the general equations are obtained from those of the diagonal case 
just replacing L by A, that can be therefore considered as a sort of 'effective 
angular momentum' in the case where the Killing fields arc not orthogonal. 

Ill Prolongability of the geodesies 

In this section we shall introduce two theorems on causal geodesic completeness 
of orthogonally transitive cylindrical spacetimes. Null coordinates, 

(19) 



t + r 
u= ^ , V- 


t-r 
2 


will play an important role in the results. 





Theorem 1: An orthogonally transitive cylindrical spacetime endowed with a 
metric whose local expression in terms of C^ metric functions /, g, A, p is 
given by (|l|) such that the axis is located at r = has complete future 
causal geodesies if the following set of conditions is fulfilled: 

1. For large values of t and increasing r, 

(a) I K>0 
[ 9u > 0, 

{5,. > or \gr\<9u 
hr > or |/v| < hu 
qr>0 or \qr\<qu- 

2. For L ^ and large values of t and decreasing r, 

(a) Sg, + P^e^fq, + A^^k,>0 

(b) Either Sgr+P^e^f qr+A^^^ /i^ < or Sgr+P^e^f g^+A^^ K < 
Sg^,+P^e^fq,+A^^K. 

3. For large values of the time coordinate t, constants a, b exist such 

2g{t,r) ] 

that g{t,r) + f{t,r)+ In p-\n\A\ \ > -\n\t + a\ + b. 
g{t,r)-f{t,r) J 

4. The limit lim — exists. 

r^O p 

A theorem can be obtained for past-pointing geodesies just changing the sign 
of the time derivatives in the previous one. 

Theorem 2: An orthogonally transitive cylindrical spacetime endowed with a 
metric whose local expression in terms of C^ metric functions /, g, A, p 
is given by (|l|) such that the axis is located at r = has complete past 
causal geodesies if the following set of conditions is fulfilled: 

1. For small values of t and increasing r, 

{ 9v<Q 
(a) I K<0 

{ qv < 0, 

{g,. > or Iffrl < -gv 
hr > or \hr\ < -K 
qr>Q or \qr\<-qv 

2. For L 7^ and small values of t and decreasing r, 

(a) Sgu + P^e^f ?„ + A^^h^ < 



p" 



(b) Either dgr+P'^e^f qr+A^^ hr<OoT Sgr+P^e^^ qr+A^^ K < 
\Sg^ + P^e^fqu + A^^K\. 



3. For small values of the time coordinate i, constants a, b exist such 

that g{t,r) + f{t,r)+lnp-\n\A\ > > - In |t + a| + 6. 
git,r)~fit,r) J 

4. The limit lim — exists. 

r^O p 

The theorems in B can be seen to be subcases of the ones introduced here. 



IV Proof of the theorems 

In order to achieve prolongability of causal geodesies we just have to impose that 
t remains finite for finite values of the affine parameter. The radial velocity, r, 
need not be considered since it cannot become singular if t is not singular too. 
The other derivatives, i and </>, are quadratures of smooth functions of t and r 
and therefore they only may turn singular if i or r become so. We shall focus on 
future-pointing geodesies. The analysis for past-pointing geodesies is entirely 
similar. 

A way of preventing the hyperbolic functions of S, from becoming singular 
is to require that ^ does not grow indefinitely. Therefore for large values of 
^, ^ must eventually become negative. Since the constants, L, P, i5 may vanish 
independently one is not to expect compensations between them. Hence their 
respective terms in (ll^) have to become negative for large values of positive ^ 
as it is stated in conditions ( |la| ) and (lb), taking into account that coshf 



sinh^ + e~^ and that therefore the terms in the negative exponential of ^ (lb) 
need not be negative but just of the same order as those in (|la|). 

By imposing condition (Q) in the theorem, we require that the geometry 
of the spacetime in the vicinity of the axis is determined by p and not by A. 
Hence for negative decreasing ^ the axis could be singular only for geodesies 
with L y^ 0, since the terms A/p are finite at the axis. Geodesies with zero 
angular momentum just cross r = and reappear with positive ^ and polar 
angle 4> + tt and hence need not be taken into account. In this case we can 
therefore allow compensations between the non-zero L term and the other ones. 
Sphtting cosh^ as e^ — sinh^, the condition for ^ to become positive for large 



values of t and negative ^ is stated in (|2a|) for the terms in sinh^ and in (2b) 



since they are exponentially damped. 

No further conditions need be imposed on ^. But t could turn singular also 
for the e~^^F term. A way of preventing it is to impose a growth slower than 
linear for t due to each of the three terms {6, A, P) for large values of t. This 
is done in condition (S). 

The condition on g must be refined since we have not yet considered the 
geodesies that cannot be parametrized by i^. These are those with zero F, 
that is, null geodesies with zero i and 0. Since i = \r\ for such geodesies, the 



equations of geodesic motion, 

i + gti^ + 2grir + gtf'' ^Q, (20) 

r + grP + 2 gtif + grV^ = 0, (21) 

can be reduced to a single one that can be integrated, 

i + 2 gtP + 2 g^ir = ^ (e^^t) ' = if, (22) 

which, in order to have t extendible to arbitrary values of the affine parame- 
ter, can be controlled by imposing at most linear growth for t as it is done in 
condition (0). 

V Completeness of several cylindrical models 

Since all known diagonal cylindrical perfect fluid models ( ||ri| , ll^ , p3| , ||I3 ) with 
regular curvature invariants have already been shown to be geodesically com- 
plete B, we shall only concern about nondiagonal ones. To our knowledge 
there are just two and both can be derived from Einstein spacetimes using the 
Wainwright-Ince-Marshman generation algorithm for stiff perfect fluids |15| . 

1. Mars: p6[ ft is the first known nonsingular nondiagonal cylindrical cos- 
mological model in the literature. In another context it was previously 
published by Letelier p^ . In isotropic coordinates the metric functions 
can be written as, 

g{t,r) = i In cosh(2 at) -I- ^aa'^r'^, fit,r) — ^ lncosh(2 a t), 

p(t,r) = r, A{t,r)^ar\ (23) 

where a is a constant and a > 1 . If a = 1 the pressure and the density of 
the fluid vanish. 

All functions are even in t and therefore the analysis of past-pointing 
geodesies is equivalent to that of future-pointing ones and shall be omitted. 

The derivatives g„, Qu in condition (|la|) are positive for positive t whereas 
hu also needs large radial coordinate. 



The derivatives gr, qr in condition (lb) are positive regardless of t and hr 
is positive for large r. 



Conditions (|2a|), (2b) are satisfied when r is small and t is positive since 
the L term is dominant for small r. 

The functions in condition (H) are all positive except for the Inr term in 
h when r decreases. But this can be bounded by a logarithm of t and 
therefore the condition is fulfilled. 



The ratio Aj p tends to zero for decreasing r and hence condition (H) in 
theorem 1 is satisfied. 

Hence this spacetime is causally g-complete. 

2. GrifHths-Bicak: fll] The previous model is comprised in this one for c = 
after a redefinition of constants. The metric functions can be written as, 



g(t,r) = ilncosh(2at) + \a^r^ + \^{t,r), f{t,r) = ^ lncosh(2 ai), 

p{t,r)^r, A{t,r)^ar^, (24) 



where J7 is a function that is obtained from a solution, a, of the wave 
equation, 

a(i,r)=M + V2c^V^^_^^;^i_^p^_^p . (25) 

The analysis of the geodesies of this spacetime can be dealt with easily 
since, 

Qu = r{at + UrY, fit, = -r(crt - o-^)^, (26) 

and therefore fJ contributes with an additional positive term to conditions 



(la), (lb) in theorem 1, that were already checked for Mars spacetime. 



On the contrary, Q. adds a negative term in conditions (Eq), (2b) but it is 
negligible for small r. 

Finally, the contribution of fi to g grows at most as y/t for large t and is 
therefore negligible compared to the other terms in condition (JSJ). 

A similar reasoning is valid for theorem 2, although the metric function fJ 
is not even in time. 

Hence these spacetimes are geodesically complete. 
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